DuulA'A ^ N 2,pi*ivDU&TE SCHOOL 

c ^ F0Rt)1A 939 



NAVAL 



POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



ONE 


DIMENSIONAL MODEL 


HINDCASTS OF WARM 


ANOMALIES IN THE NORTH PACIFIC OCEAN 




by 






Bernardino J. 


Jaramillo 




December 


1984 


Thesis 


Advisor : 


R.W. Garwood 



Approved for public release; distribution unlimited 




UNCLASSIFIED 

SECURITY CLASSIFICATION of THIS PAGE (When Dmtm Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subtitle) 

One Dimensional Model Hindcasts of Warm 
Anomalies in the North Pacific Ocean 


5. TYPE OF REPORT 4 PERIOD COVERED 

Master's Thesis 
December 1 984 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHORfaJ 

Bernardino J. Jaramillo 


8. CONTRACT OR GRANT NUMBERf«J 


9. PERFORMING ORGANI ZATION NAME AND ADDRESS 

Naval Postgraduate School 
1 Monterey, CA 93943 


10. PROGRAM ELEMENT, PROJECT TASK 
AREA 4 WORK UNIT NUMBERS 

1 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, CA 93943 


12. REPORT DATE j 

December 1984 


13. NUMBER OF PAGES I 

8 2 j 


1 14 MONITORING AGENCY name & A D D R ESSf If different from Controlling Office) 


15 SECURITY CLASS (ol thta report) 

Unclassi f ied 


1 5*. DECLASSIFICATION/DOWNGRADING 

SCHEDULE | 



16. DISTRIBUTION STATEMENT (ol this Report) 



Approved for public release; distirubit ion unlimited 



17. DISTRIBUTION STATEMENT (of the abstract entered In Block 20, If different from Report ) 



18. supplementary notes 



19. KEY WORDS (Continue on reverse side If necea sary and identify by block number) 

North Pacific Ocean 
Warm Anomalies 

One-Dimensional Model Hindcasts 

TRANSPAC 

20. ABSTRACT (Continue on reverse side If necesaary and Identify by block number) 

Hindcasts of the development, maintenance and decay of warm 
anomalies in the North Pacific Ocean are attempted with the 
Garwood one-dimensional ocean mixed layer model. Monthly 
temperature analyses from the North Pacific Ocean Experiment 
(NORPAX) provide the initial and verifying fields, and the 
atmospheric forcing is derived from FNOC archived analyses. 



DD 



FORM 
1 JAN 73 



1473 



EDITION OF 1 NOV 65 IS OBSOLETE 

S N 0102- LF- 014- 660 1 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS PAGE (When Dmtm Sntmrmd) 



UNCLASSIFIED - 



SECURITY CLASSIFICATION OF THIS PAGE flTh#n D«<« En<»r» <S) 



Large systematic errors in the anomaly hindcasts during the early 
summer are apparently due to erroneous surface heat flux 
specifications. Variations of the heat flux corrections proposed 
by Elsberry et al. (1982) are extensively tested for one warm 
anomaly development period. The climatology of Wyrtki (1967) is 
the heat flux specification most successful in simulating the 
observed anomaly patterns. The hindcasts demonstrate 



importance of 
prediction of 



an accurate heat 
warm anomal ies . 



flux specification for 



the 

correc t 



S ' N 0102- LF- 014- 6601 



2s 



UNCLASSIFIED 



ECUR1TY CLASSIFICATION OF THIS P AG EpITian Dmtm Entmrmd) 



Approved for public release; 



distribution is unlimited. 



One Dimensional Mcdel Hindcasts 
of Warm Anomalies 
in the North Pacific Ocean 



by 



Bernardino J. Jaramillo 
Lieutenant Commander, United States Navy 
B.S., University of Utah, 1975 



Submitted in partial fulfillment of the 
requirements xor the degree of 



MASTER OF SCIENCE IN METEOROLOGY AND OCEANOGRAPHY 



from the 



NAVAL POSTGRADUATE SCHOOL 
December 1984 



ABSTRACT 







Hindcasts of the development, maintenance and decay of 
warm anomalies in the North Pacific Ocean are attempted with 
the Garwcod one- di mensional ocean mixed layer model. Monthly 
temperature analyses from the Ncrth Pacific Ocean Experiment 
(NORPAX) provide the initial and verifying fields, and the 
atmospheric forcing is derived from FNOC archived analyses, 
large systematic errors in the anomaly hindcasts during the 
early summer are apparently due to erroneous surface heat 
flux specifications. Variations of the heat flux corrections 
proposed by Elsberry et al. (1982) are extensively tested 
for one warm anomaly development period. The climatology of 
Wyrtki (1967) is the heat flux specification most successful 
in simulating the observed anomaly patterns. The hindcasts 
demonstrate the importance of an accurate heat flux specifi- 
cation for correct prediction of warm anomalies. 
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I. INTRODUCTION 



A. PURPOSE AND HYPOTHESIS 

The purpose of this study is to test the applicability 
of the Garwood (1977) one-dimensional bulk mixed layer model 
for hindcasting the warm ocean thermal anomalies found in 
the Anamoly Dynamics Study (ADS) domain. Descriptions of 
several warm anomalies suitable for this study are found in 
Elsberr y (1 984) . 

Stringer (1983) conducted a similar study which involved 
cold anomalies in the ADS domain. He found that the Garwood 
model does a good job of hindcasting cold anomalies in the 
fall-winter period but experiences larger errors in the 
spring-summer period. Stringer suggests that uncertainties 
in the heat flux correction fields are fairly large in the 
spring and summer and limit the ability of the model to 
hindcast anomalous ocean conditions during those seasons. 
Therefore, a second purpose for this study is to further 
examine the heat flux correction fields and to uncover 
factors which may ' affect the ability of the model to hind- 
cast in the spring and summer. 

Three warm anomalies have been chosen for this study 
that should satisfactorily test the Garwood model for the 
spring and summer conditions. The first case involves a 
warm anomaly starting in January 1976 and ending in July 
1976. The next two anomalies extend from June 1977 to 
December 1977 and June 1973 tc November 1978. The time 

periods of these anomalies should adequately test all 
aspects of the Garwood model for the spring and summer 
hindcasts. 
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B. STUDI DESCRIPTION 

The North Pacific Ocean Experiment (NORPAX) ADS area 
(Figure 1.1) is the region examined in this study. The 
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Figure 1.1 ADS REGION. 

region is bounded by 30N-50N and 130W-160E. Analyses of 
ships-of-opport unity expendable bathythermographs (XBT) 

provide ocean thermal structure fields on space scales of 
thousands of kilometers and time scales of one month (White 
and Bernstein, 1979). The TRANSPAC analyses provide the 
initial data for the ocean prediction model and also serve 
as validation fields. 

As briefly described in section A, three warm anomaly 
periods are examined in this study. The first set of warm 
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anomalies (WA 1/2) occurs during winter (January 1976) and 
continues through July 1976. The second set of warm anoma- 
lies (WA 12/13/15) are part of a large-scale pattern which 
developed during the summer (June 1977 to December 1977). 
The third set from June 1978 to November 1978 contains an 
extreme, near-surface warm anomaly (WA 23) that forms near 
the western boundary of the domain. All three anomaly cases 
provide excellent examples for the model hindcast studies 
recommended by Stringer (1983). 



II. DATA sources and preparation 



A. DATA sources 



Fleet Numerical Oceanography Center ( 
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B. PBEDICTION MODEL 



The model used for this study is the Garwood (1977) 
second-order closure, bulk model. Advection is not included 
in the model. The Garwood model is initialized with THANSPAC 
objectively analyzed temperature fields, which are assumed 
to apply on the 15th of the month. Model output consists of 
vertical temperature profile predictions in 5m increments to 
200m over the entire grid every hour. These hourly tempera- 
ture profiles are averaged over five days centered on the 
15th of each month. An additional output is the mixed layer 
temperature at the time of maximum mixed layer depth for 
each day. Model anomaly fields are formed by subtracting the 
four-year mean (climatological) fields from the predicted 
fields at 0, 60 and 120m. 



III. HASH ANOMALY DESCRIPTION 



A. HARM ANOMALY 1/2 

The warm anomalies (WA) described in this section were 
initially designated and described by Elsberry (1984) as WA 1 
and WA2. This case persisted for approximately six months, 
as shown in Figures 3.1 through 3.6. 

W A 1 is evident to a depth of 60m during the first month 
of the period (Figure 3.1) in the southeast portion of the 
region. The warm anomaly intensifies to maximum amplitude 
(2.02 C at 34N, 1 5 0W) at the surface in February (Figure 

3.2) . It decreases in intensity as the period progresses 
(Figure 3.3 and 3.4). In May, W A 1 (Figure 3.5) disappears at 
the surface, but it is still evident at 60m. In June(Figure 
3.6) , a warm anomaly again appears at the surface in rela- 
tively the same position, but it is unclear whether this is 
still part of WA 1 or is an entirely new feature. 

WA2 is barely evident at approximately 38N,170E during 
January (Figure 3.1). It continues to intensify during the 
period. In March (Figure 3.3) , it has become a significant 
warm anomaly extending to 120m. In April and May (Figures 
3.4 and 3.5), WA2 appears to Hove north and elongate while 
losing some of its structure at depth. Toward the end of the 
period (Figure 3.6), WA2 becomes disorganized at the surface 
but maintains some structure at depths of 60 and 120m. 



B. -HARM ANOMALY 12/13/15 



Warm anomalies 12/13/15 were initially designated and 
described by Elsberry (1984) . This case continued for 
approximately six months (Figures 3.7 through 3.12). 



WA 12 develops in July 1977 in the eastern part 
(38N,150W) of the region and reaches an anomaly magnitude of 
+1.7C (Figure 3.7). Development continues in August (Figure 
3.8) and appears to have reached a depth of approximately 
120m, although the anomaly center is shifted slightly east- 
ward at 60 and 120m. During September and October, FA 12 
weakens and disappears (Figures 3.9 and 3.10). 

WA 13 also develops in July 1977 in the region around 
48N, 170W (Figure 3.7). In August 1977, its position and 

movement appear to be influenced by a large cold anomaly at 
approximately 40N, 170W. This ccld anomaly seems to cause WA 

13 to merge with another warm anomaly to the southeast of 
the initial position of WA 13 (Figure 3.8). Evidence of the 
merging in September (Figure 3.9) is indicated by the large 
broad warm anomaly in the northeast portion of the region. 
The large anomaly diminishes in size in October and 
November, and it is no longer evident in December (Figures 
3.10 and 3.11). 

WA 15 is evident in July 1977 (Figure 3.7) along the 
extreme western boundary at 40N, 160E. It continues to 
develop and reaches maximum imtensity of + 3.63C in August at 
42N, 170E (Figure 3.8). At this time it also appears to have 

merged with WA 13, which contributes to a large warm anomaly 
area along the northern border. The anomalously high temper- 
atures persist through November and disappear in December 
(Figures 3.11 and 3.12). 

C. WARM ANOMALY 23 

Warm anomaly 23 was initially designated and described 
by Elsberry (1984). This case persisted for approximately 
four months as seen in Figures 3.13 through 3.16. 
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During June 1978 (Figure 3. 13), a general warm area at 
the surface has developed on the southwestern boundary of 
the region. By July (Figure 3. 14), this area shows a strong 
(approximately +2.5C) surface anomaly at the boundary 
region. A slight shift northward is evident in August 
(Figure 3.15), and the temperature anomaly decreases by 
approximately 1.5C. By September (Figure 3.16),WA 23 has 

disappeared . 
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Figure 3. 1 Temperature (C) anomaly WA 1/2 during 
January 1976 at (A) surface, (B) 60m and (C) 120m. 
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Figure 3.2 Similar to Figure 3. 1 except for February 1976. 
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Figure 3.3 Similar to Figure 3.1 except for March 1976. 
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Figure 3.4 Similar to Figure 3.1 except for April 1976. 
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Figure 3.5 Similar to Figure 3.1 except for May 1976. 




Figure 3.6 Similar to Figure 3.1 except for June 1976. 
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Figure 3.7 Temperature (C) anomaly HA 12-15 for July 
1977 at (A) surface (B) 60m and (C) 120m. 
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Figure 3.8 Similar to Figure 3.7 except for August 1977. 
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Figure 3-9 Similar to Figure 3.7 except for September 1977. 
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Figure 3.10 Siailar to Figure 3.7 except for October 1977. 
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Figure 3.11 Similar to Figure 3.7 except for November 1977. 
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Figure 3.12 Similar to Figure 3.7 except for December 1977. 
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Figure 3. 13 Temperature 
1978 at { A J surface 



(C) anomaly HA 23 
(B) 60m and (C) 
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120m. 
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June 
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Figure 3.14 Similar to Figure 3.13 except for July 1978. 
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Figure 3.15 Similar to Figure 3.13 except for August 1978. 
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Figure 3.16 Similar to Figure 3.13 except for September 1978. 
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IV. MODEL HINDCASTS 



A. INITIAL HINDCASTS 

Several model integrations were cond 
hindcast ability of the Garwood model for 
cases described earlier. The correction 
the FNOC forcing fields were described 
(1982). The correction fields, which 
assuming a local heat balance during a 
consist of six bi-monthly fields which a 
the FNOC heat fluxes. The corrections are 
heat fluxes to be used for one- dimensional 
These corrections eliminate a bias of 
surface heat flux found in the FNOC fields 



1 . Model Hindca sts of WA 1 /2 
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Comparisons with Figures 3.4 through 3.7 show that there is 
good agreement in April at a depth of 60m, but the placement 
of the features at the surface and 200 m does not agree with 
the observations. Similar comparisons for May and June show 
even less agreement as the model continues to warm exces- 
sively. 

An additional model hinacast was conducted using May 
as the initial month. Again, the hindcast was too warm, and 
the results very similar to the June initialization. To 
determine the effect of the correction to the forcing field, 
two model integrations were conducted from January through 
July and May through August without application of the 
correction field. Neither of these model hindcasts agreed 
with the observed anomaly fields in that both predicted a 
warm anomaly in the northwest area and a large area of 
anomalously low temperatures in the entire southern region 
from approximately 40N southward. Further studies with modi- 
fied correction fields in space and time will be described 
in section C of this chapter. 

2. Model Hindcasts of WA 12/13/15 

Model hindcasts of WA 12/13/15 were generally cooler 
than observations over much of the area with a pronounced 
warm region in the southwest portion of the region. Using 
the June analysis to specify the initial conditions of the 
model (Figure 4.5), it is seen that July (Figure 4.6 ) is 
already much cooler than climatology in the northeast and 
much warmer in the southwest. Comparison with observations 
for July (Figure 3.7) shows very little correlation between 
the two plots. The August plot (Figure 4.7) maintains the 
same pattern developed in July and does not bear any resemb- 
lance to the observed July condition (Figure 3.8). 
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The model was then tested with initial conditions 
from March and May to determine if a relationship existed 
between initialization period and accuracy of model hind- 
casts. In both cases, the model predicted development of a 
general warm area in June, although by July there were cold 
anomalies in the northeastern portion of the area. 3y 
August, the size of the warm area was increasing (Figures 
4.8 and 4.9). This was similar to the model output with a 
June initialization. Thus, time of initialization does not 
appear to affect the output as there are anomalously low 
temperatures as early as July. The model was also integrated 
without a correction to the forcing field. The output (not 
shown) for a model hindcast from June inital conditions 
showed the entire grid area was anomalously cool in July, 
and the temperatures continued to decrease as the period 
progressed. The different variations used to hindcast WA 
12/13/15 did not significantly improve any of the hindcasts. 
Thus, the early summer period appears to pose serious prob- 
lems for reliable model hindcasts. 



3. Model Hindcasts of WA 2 3 



Two model integrations 
initialization time of January 
zation of June (Figures 4.10 th 
casts produced a large anomalo 
the region with a small warm an 
As the period progressed throug 
cold anomalies at the surface 
unverified cold pattern. C 
(Figures 3.13 through 3.16) in 
feature placement or feature 
Temperature structure with 
favorably. 



for WA 23 were done with an 
and June. With an initiali- 
rough 4.12) , the model hind- 
usly cool area over most of 
cmaly on the southern border, 
h August and September, the 
intensified and developed an 
cmparison with observations 
dicate a lack of agreement in 
intensity during all months, 
depth also did not agree 
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B. DISCUSSION 



General results of each model hindcast showed that the 
model anomalies for VIA 1/2 were hindcast well until April. 
However, the model tended to produce an anomalously warm 
condition throughout the remainder of the period. Model 
anomalies for V? A 12/13/15 were hindcast as too cold as early 
as the second month of the model integration. Model anoma- 
lies for VIA 23 were also hindcast as too cold throughout the 
period. 

The three cases have at least one common factor in that 
they span the early summer months of June and July. During 
these months, changes in forcing (heat fluxes and wind 
speeds) are much more difficult to predict and can not be 
corrected by a bi-monthly average correction technique as 
previously done (Stringer, 1983) for cold anomaly studies. 
The variability in the total heat flux can cause the model 
to over-correct or under-correct in the quantity of heat 
necessary to simulate accurately the observed conditions for 
the early summer. 

The importance of correctly determining the heat storage 
and heat input is discussed by Wyrtki and Ulrich ( 1 382) . 
They found that over a period of one month the observed 
change in heat storage can be computed to only an accuracy 
of +/- 100 W/m 2 , which is as large as the net daily heat 

exchange. This uncertainty makes it difficult to account 
for contributions of advection or mixing over a monthly 
interval. If this magnitude cf uncertainty is encountered 
in actual ocean estimates, then it is very possible that 
errors this large may exist in the FNOC short-term 
predictions used here for the forcing fields. Uncertainty 
in determining the heat storage and heat input can also 
adversely affect the depth of the mixed layer as predicted 
by the model. For a shallowing mixed layer depth with no 
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entrainment, the mixed layer depth , 11 , is related to the 
total heat flux, Q c , by the relationship: 

h-Ly u* 3 / C^gQo / (fC r ) } where 
L=Obukhov length scale 
u*=Friction velocity 
Cp=Specific heat of water 
Q c =Heat flux 
<*■= Specific volume 
f =Density 
g=Gr avity 



The corresponding temperature change is given by 

dT/dt= (Q„//g)/h. 

The heat flux may be expressed as 

Qo -Qo (correct) +Qo (error) =Qc +Q e =Q C (1 +E) , where 
E is an error bias equal to Q e /Q c . This now gives 

h=fC p u*3 / {^gQc ( 1 + 2) } = h c / (1 + E) 
and with the expression for Q c yields 

dT/dt= 2 (i +E ) ^g/{ (/9) 2 U *3 } 
d I/dt= { (1+E) 2} {dT/dt} c 

This expression demonstrates that a net flux error, Q , 
contributes as a square because it also contributes to an 
incorrect mixed layer depth. 

The heat flux corrections by Eisberry et al. (1S82) were 
computed on a bi-monthly basis and were averaged over three 
years. This average correction field does not appear to 
accurately modify the forcing field in the proper amounts 
for the early summer. To study this problem further, the 
correction field for WA 1/2 was modified, and a different 
source of heat flux values was used to improve model 
hindcasts. 
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C. TJ A 1/2 DETAILED STUDY 



1 . One -Dimensio nal Heat Budget Correction Field 



Previous case studies of cold anomalies by Stringer 
(1983) and the initial model hindcasts for warm anomalies in 
this study indicated that the model hindcasts are degraded 
by boundary condition (surface heat flux) errors in the late 
spring and early summer. Since the bi-monthly average 
correction field derived by the method of Slsberry et al. 
(1982) did not correctly modify the forcing fields for late 
spring and early summer cases, a new correction field was 
derived. This correction field exactly accounts for the 
difference between the heat content in the water column to a 
depth of 200 m and the integrated surface heat at a given 
geographical location. The correction is formed from the 
following: 



I=jQr(t) at 
h o 

H= /6VJ(T (z) -T (200m) } dz 



I is the integrated surface heat flux, H is the heat 
content, Q r is the surface heat flux, t is time, t) is the 15th 
of one month and tjis the 15th cf the following month. This 
yields the unfiltered exact correction (<ff) ; 

£ = {H ft) -H ty } -I*/ ( h ' t ! ) h curs 

This field was then filtered using the same procedureas in 
the original correction field. 

The one-dimensional heat budget correction field was 
used in nodel hindcasts for May through August 1976 . These 
hindcasts were compared with monthly temperature change 
fields. This comparison emphasizes model temperature 
changes and is not dependent on climatology as in the 
anomaly comparisons above. 
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For June to May model temperature fields (Figure 
4.13) are close to the observed fields (Figure 4.16). In 
overall structure, the field shows some similarity to the 
observed pattern, especially in the warming of approximately 
2C in the eastern third of the region. 

The July to June model plot (Figure 4.14) shows two 
distinct warm features at 42N,175S and 32N,165W which 
deviate greatly from the observed fields (Figure 4.17). 
Larger temperature difference gradients are more evident in 
the model plot. Locations of maximum warming and cooling do 
not correspond to similar features in the observed field, 
and it appears that the model is generating temperatures too 
high for the entire region during July. 

The August to July (Figure 4.15) model output 
appears to have even less correlation with the observed 
field (Figure 4. 18) than during the previous two months. The 
two general areas with large deviations noted previously 
again appear. The overall pattern of large gradients in 
both warming and cooling constitutes an unverified hindcast. 

To determine what caused the model to deviate so 
greatly from observations in July and August, observed and 
model temperature profiles at 42M, 175E were analysed 
(Figures 4.19 and 4.20 ). It is seen from the June and July 
observed temperature profiles that a large amount of heat is 
added in the water column during this two-month period. 
Calculations of the heat storage show that approximately 
3.0x10 4 cal/cm 2 is added, which is abnormally large when 
compared to surrounding values. Analysis of the surface heat 
flux field shows a value of I.OxlO 4 cal/cm 2 which is in good 
agreement with surrounding points.- Therefore, it appears 
that the large heat content icfluenced the correction and 
caused the model to over-correct by adding too much heat. 
The reason for such a large heat content could be horizontal 
advection of warmer water, downward vertical velocities, or 



39 



possible observational error. Whatever caused this heat 
excess is not handled well by the model, and must be better 
understood to be corrected. Use of the or.e-diraensional heat 
budget correction field is not sufficient for the model to 
produce acceptable hindcasts for all points. 

2. Correction to the One- Dimensional Heat B udge t 

Correcti on F iel d 

In an attempt to correct the large deviations seen 
in the model using the exact correction field, a correction 
method applied to a specific area was devised. The point 
(42N,175E) was chosen for this test since there is a large 
temperature difference between July and June and the model 
temperature profile is greatly exaggerated for July. The 
point correction, £ P , was derived by removing the amount of 
heat in the mixed layer during the two-month period and 
using the residual heat content as an additional correction 
factor to the correction field at that point. For (42N,175E) 
the following calculations were used: 

£p =AH ZCC -AH (7 > l /£> =28. 1 3 cai/cm 2 -hr 

where A H zcc =3x10 4 cal/cm 2 and A = 1.0x10 3 
cal/cm 2 and £ = 720 hours. 

This correction was then applied to the one- 

dimensional heat budget correction field at that point, and 
the model was integrated again. The resulting model tempera- 
ture profiles (Figure 4.21) are much closer to the observed 
temperature profiles. The July profile is within 0.5C at the 
surface and is in close agreement down to 80 m. The August 
model profile still shows a large amount of heat in the 
column and requires some additional correction to reproduce 
the observed temperature profile. 

Although this procedure did improve the July temper- 
ature profile, it would not be useful to apply a similar 
correction at every point without a better understanding of 
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the cause of the error (i.e., neglect of advection or obser- 
vational error) . A better method for deriving either the 
correction field or the initial forcing field is needed to 
make better use of the model as a prediction tool. 



3 . Average Cl im ato logic al Surface Heat Flux 
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the solar radiation 



was added as a separate field that was also derived from the 
climatological charts of Wyrtki. Thus, below-surf ace 
absorption was considered in Case II. Case III was similar 
to Case II, except it included the diurnal cycle for the 
radiation. 

C as e I: Plots of monthly temperature changes for 

June to May (Figure 4.30) were much smoother than the 
previous results using the' FNOC heat fluxes (Figure 4.13) 
and temperature magnitudes were much closer to the observed 
values (Figure 4.16). The region has a general warming 
pattern increasing from west to east and north to south with 
a small temperature change maxinum in the southeast corner. 
The July to June plot (Figure 4.31) shows a larger gradient 
from north to south with a much more intense warming area 
situated at approximately 32N,170W. The pattern tends to be 
warmer than observations (Figure 4.17) with the largest 
temperature changes shifted slightly south and east of a 
similar but lower temperature change maximum in the observa- 
tion. The August to July plot (Figure 4.32) also has a 
smoother field and shows a slight decrease in temperature 
change. The bull’s-eye pattern apparent in the previous plot 
has been modified and appears as a broad warm area in the 
southwestern part of the grid. These temperatures are 
slightly lower than observations (Figure 4.18) . General 
warming throughout the entire region is still evident. 

Although placement of features is not completely 
consistent with observations, the general overall warming 
pattern is verified. Temperature change magnitudes are 
close to observed values. Temperature patterns produced by 
the model tend to be somewhat zcnal due to the zonal pattern 
of the climatological heat flux fields. 

Case II: Inclusion of a separate solar radiation 

field did modify the temperatures compared to Case I. The 
June to May plot (Figure4. 33) is very similar to that in 
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Case I with a reduction of the temperature m 
approximately 1.0C. The July to June plot (F 
shows a similar trend, with temperatures reduc 
approximately 1.0C. The August to July plot ( 
also shows a pattern similar to that seen in 
slightly smaller temperature changes. 

Overall, the inclusion cf a separate sol 
field produced patterns very similar to that in 
gradients of temperature difference were much sm 
the observed general warming trend is correctly 
Case III: The diurnal variation to the 

tion field further modified the temperature 
pattern seen in Case II. All three plots (F 
through 4.38) are very similar in structure to 
plots previously discussed. The July to June pi 
greatest modification with an additional decreas 
the area of the bull’s-eye described earlier, 
warming pattern is still apparent and, even with 
cation of temperatures and weaker gradients, 
similar to that of the observations is predicted 
Climatological heat flux fields do aid 
producing patterns which are closer to observati 
previous attempts using FNCC heat flux fields 
placement of features seen in the observations 
handled well by the use of climatology for th 
forcing . 
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Figure 4.1 Model teaperature anomaly MA 1/2 during April 1976 
at (A) surface, (B) 60a, and (C) 120m- 
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Figure 4.2 Siailar to Figure 4.1 except for May 1976. 
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Figure 4.3 Similar to Figure 4.1 except for June 1976. 
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Figure 4.4 Siailar to Figure 4.1 except for July 1976. 
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Figure 4.5 Model temperature anomaly MA 12-15 during 
June 1977 at (A) surface, (B) 60m, and (C) 120m. 
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Figure 4.6 Similar to Figure 4.5 except for July 1977. 
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Figure 4.7 Similar to Figure 4.5 except for August 1977. 
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Figure 4-8 Siailar to Figure 4.5 except for 
model initialization in May 1977. 
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Figure 4.9 Similar to Figure 4.8 except for July 1977- 
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Figure 4.10 Model temperature anomaly 
at (A) surface, (B) 60m, and 



MA 23 during July 1978 
(C) 120m. 
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Figure 4.11 Similar to Figure 4.10 except for August 



1978 
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Figure 4.12 Similar to Figure 4.10 except for September 1978. 



55 





Figure 4.13 Model Monthly temperature difference MA 1/2 
for June to May with the one-dimensional 
heat budget correction field. 
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Figure 4.14 Similar to Figure 4.13 except for July to Jun 
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Figure 4.15 Siailar to Figure 4.13 except for August to July, 
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Figure 4.16 Observed Monthly temperature difference 

HA 1/2 for June to Hay. 
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Figure 4-17 Siailar to Figure 4.16 except for July to June- 
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Figure 4.18 Similar to Figure 4.16 except for August to July. 
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Figure 4.19 Hodel temperature profiles for May through August 

1976 at 42N, 175E with 
one- dimensional heat fcudget correction. 
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Figure 4.20 Similar to Figure 4.19 except for observed values. 
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Figure 4.21 Similar to Figure 4.19 except with added corrections. 



64 



SON 
45N 
40N 
35N 
30N 

1 60E 1 70E 1 80 1 70W 1 60 W 1 SOW MON 1 30N 




Figure 4.22 Climatological total heat flux for May. 




Figure 4.23 



FHOC average total heat flux for May. 




Figure 4.24 Siailar to Figure 4.22 except for June. 







Figure 4.25 Siailar to Figure 4.23 except for June. 
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Figure 4.27 Similar to Figure 4.23 except for July. 
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Figure 4.29 Similar to Figure 4.23 except for August. 
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Figure 4.30 Model monthly temperature difference 
for June to May with climatological heat flux field 

69 



MA 1/2 
(Case I) . 




Figure 4.31 Similar to Figure 4.30 except for July to June. 
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Figure 4.32 Similar to Figure 4.30 except for August to July. 
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Figure 4.33 Similar to Figure 4.30 except with a 
climatological radiation field (Case II) . 
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Figure 4.34 Similar to Figure 4.33 except for July to June. 
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Figure 4.35 Similar to Figure 4.33 except for August to July. 
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Figure 4.36 Similar to Figure 4.30 except with a 
diurnal radiation field (Case III) . 



75 




-I - r - “ - i 1 1 1 1 

1 6 0 E 1 70E 1 80 1 7 0 W 160W 150W 1 40W 130W 



ure 4.37 Similar to Figure 4.36 except for July to June. 
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Figure 4.38 Similar to Figure 4.36 except for August to July. 



77 



V. SUMMARY AND CONCLUSIONS 



None of the three warm ano 
this study were hindcast well b 
model using heat fluxes derived 
hindcasts seemed to develop d 
summer, one case (WA 1/2) in 
study . 

A detailed examination of 
fields used in previous studies 
1983) revealed that the early s 
a seasonally varying heat flux 
the average heat flux fields di 
warm anomalies as for the earl 
local heat budget correction 
improve hindcasts. Some impr 
areas was seen when selective 
this procedure would be of 
predictions. For research pu 
correction is a useful method 
difficulties in the hindcast mo 

Use of climatological heat 
with FNOC winds produced a bett 
FNOC's heat fluxes. FNOC hea 
various forms used in this stu 
which more closely resembled ob 
ture anomalies are primarily at 
mixing, the use of climatolog 
winds (synoptically varying) 
casts of thermal structure. Ho 
used in this study, still doe 
layer dynamics attributable to 



maly (WA) periods chosen for 
y the Garwood one-dimensional 
frnm FNOC. Since all three 
ifficulties during the early 
1976 was selected for further 

the FNOC surface heat flux 
of cold anomalies (Stringer, 
ummer can not be corrected by 
correction. Corrections to 
d not prove as successful for 
ier cold anomaly studies. A 
field did not significantly 
ovement in specific problem 
corrections were done, but 
little value in real-time 
rposes, however, selective 
to diagnose potential model 
de of operation. 

flux fields in conjunction 
er result than did the use of 
t fluxes and winds and the 
dv produced general patterns 
servations. Whenever tempera- 
tributable to anomalous wind- 
ical heat fluxes with actual 
ray provide reasonable hind- 
wever, this climatology, as 
s not account for the mixed 
synoptic variability in the 
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heat fluxes. Horizontal and vertical aavection are also 
neglected, and may account for some of the differences 

between the observations and the model simulations. Also, 

the one month resolution afforded by the observations is not 
adequate to verify a model which is most sensitive to 

synoptic variability in the atmospheric forcing. 

This study has shown that it is possible to hindcast 
warm temperature anomalies for the mid-latitudes. The impor- 
tance of using a reasonably correct heat flux field was 
shown dramatically. This study indicates a need for much 
additional research to improve the surface heat fluxes. 
Only then will the model be able to predict (hindcast) 

.specific features that are related to the surface thermody- 
namic boundary conditions. Future three-dimensional 

modelling will include the advective processes neglected 
here, but the largest difficulty and potential source of 
error still appears to be the specifica tion of the surface 
heat flux. 
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